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Abstract. The numerical investigation of wave propagation in the asymptotic domain of Kerr 
spacetime has only recently been possible thanks to the construction of suitable hyperboloidal 
coordinates. The asymptotics revealed a puzzle in the decay rates of scalar fields: the late-time 
rates seemed to depend on whether finite distance observers are in the strong field domain 
or far away from the rotating black hole. In this paper we study late-time decay rates using 
a horizon-penetrating, hyperboloidal slicing with transmitting layers attached to a compact 
domain in Boyer-Lindquist coordinates. The technical construction of transmitting layers for 
Kerr spacetime should be useful in future studies of wave propagation. We discuss splitting 
of local decay rates in certain projected modes and in the full field. The splitting in the full 
field rates is explained by competition between projected modes. For the splitting in certain 
projected modes we argue that, asymptotically in time, the strong field rates are valid at all 
finite distances, but at any given late time, there are three domains with different local decay 
rates whose boundaries move during evolution. 



1. Introduction 

Generic wave propagation in a black hole spacetime consists of three stages: an initial 
transient during which the evolution depends on details of initial data, exponentially decaying 
oscillations called quasinormal mode ringing, and a polynomial late-time decay. In this paper, 
we focus on the late-time behavior of the evolution during which the field decays as t™, where t 
denotes time and n < denotes the decay rate (U). We resolve a recent puzzling observation 
regarding the decay rates of scalar fields in the asymptotic domain of a Kerr spacetime. 

The late-time decay of scalar fields in Kerr spacetime in the strong field domain has 
been a subject of controversy (||2] [3] |4] |5] |6l 0), which has been explained by differences 
in initial data (J8j 0). In a recent development, the decay rates given by a complicated 
formula involving the azimuthal mode number of the initial perturbation was recast by Burko 
and Khanna in a way that does not include the azimuthal mode number ( ifTOl ). Today, our 
understanding of late-time Kerr decay rates includes as an essential element the distinction 
between initial data modes and projected modes generated by coupling. Denoting the initial 
mode £' and the projected mode £, the strong field rates of projected modes proposed in 
Ref. HOI read n = -(£' + £ + 1) for £ < I' and n = -(£' + £ + 3) for £ > f . 

Until recently, the numerical study of Kerr tails in the literature has exclusively focused 
on the rates at finite distances, even though null infinity rates are more interesting for observers 
at astronomical distances ( ifTTl \Y2} ). Null infinity decay rates in Schwarzschild spacetime 
have been studied with the characteristic method as early as 1994 (|13|). The characteristic 
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method, however, is difficult to extend to Kerr spacetime (" lff4l [151 [l6l ). It is due to this 
technical obstacle that Kerr tails in the asymptotic domain could not be studied numerically. 

This problem has been resolved with the general construction of hyperboloidal, 
coordinates for asymptotically flat black hole spacetimes ([17]). These coordinates provide 
numerical access to the asymptotic domain, avoid the outer boundary problem, and lead to 
efficient computations because they are adapted to the propagation of outgoing waves. The 
first numerical calculation of Kerr tails along null infinity has been performed in Ref. |18|. 
Since then, more detailed numerical results using hyperboloidal compactification have been 
obtained in Refs. 1 19 20 1, which confirmed the simple decay formula proposed in Ref. iflOll 
and extended the analysis to the asymptotic domain. The decay rates of projected modes at 
null infinity are n = £ + 2 for i > £' and n = £' for £ < £' - 2, as predicted by Hod 0. 

The hyperboloidal evolutions revealed a puzzling behavior for the decay rates in the 
asymptotic spatial domain. The rates for certain observers at finite distances far away from the 
black hole deviated from both the horizon rates and the null infinity rates. This observation 
requires some explanation. In Schwarzschild spacetime, the local (in time) decay rates of 
finite distance observers depend on the distance of the observer to the black hole. Take the 
decay of the dominant £' = £ = mode. Along null infinity the field decays with a power of 
—2, and at finite distances with —3, asymptotically in time. One expects that at finite distances 
the local rates would vary monotonously between —2 and —3 depending on how far away the 
observer is to the black hole, and eventually approach —3. This expectation is in accordance 
with theoretical predictions (see Eq. (|7]) and [21, 22 1) and is confirmed by numerical studies 

dHUIH). 

A similar distance-dependence for the full field is observed in Kerr spacetime for 
£' = 0,1,2,3, but for the data with £' = 4 an anomalous behavior appears. Instead of 
a monotonous transition between the null infinity rate —4 and the finite distance rate —5, 
there are far away observers for which the decay rate seems to approach —7. In Ref. |18|, 
where this behavior was first observed, it was speculated that this is an intermediate property 
due to validity of Schwarzschild decay rates for the lowest excited modes far away from the 
Kerr black hole. A more detailed study of tail decay rates by Racz and Toth in Ref. |fl9l 
provided further evidence for this anomalous behavior. In addition, they discovered that 
certain projected modes also have different local decay rates close to the horizon and far 
away from it, which they referred to as splitting. Numerical computations in 3D by Jasiulek 
showed that splitting is robust (|20|). 

In this paper, we study the two types of splitting, namely splitting in the projected modes 
( |fl9l l20l ) and splitting in the full field ( lfl8l ). We argue that, asymptotically in time, the 
strong field decay rates are valid at all finite distances and splitting refers to an intermediate 
behavior of local decay rates only. At any given late time, however, there are three domains 
with different local decay rates for certain projected modes with £' = I. The boundaries 
between these domains move to larger distances during evolution. 

A technical development reported in this paper is the application of the hyperboloidal 
layer method from Ref. ||231 to the construction of horizon-penetrating, hyperboloidal 
coordinates based on a Boyer-Lindquist slicing. The causal behavior of the slicing is similar 
to the slicings of Refs. [ 19 20| with the difference that standard Boyer-Lindquist coordinates 
can be used in a compact domain. Such coordinates should be useful in future numerical 
studies of test fields in Kerr spacetime. 
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2. Methods 

2.1. Transmitting layers in Kerr spacetime 

We need to use hyperboloidal coordinates to investigate the behavior of fields in the 
asymptotic domain of Kerr spacetime. Most numerical codes, however, employ a Cauchy 
slicing for the simulations. From a practical point of view it seems unnecessary to modify the 
coordinates everywhere just to have access to the asymptotic domain. In this section, we use 
the hyperboloidal layer method fll23lD to construct suitable coordinates that leave a compact 
domain in Boyer-Lindquist form. Our coordinates are similar to those in Ref. E4| . but we 
extend them to become horizon-penetrating using a transmitting layer as in Ref. 11251 . 

The hyperboloidal layer method of Ref. ||231 is analogous to absorbing (or perfectly 
matched) layers in applied mathematics ([26 1) in that it attaches a layer to a finite 
computational domain in which different coordinates are used for the numerical computation. 
Instead of absorbing the outgoing waves, however, the hyperboloidal layer transmits them 
to null infinity mapped to the outer boundary of the layer (hence it is also referred to as a 
transmitting layer). There are no reflections from the layer in the continuum limit because the 
interface between the layer and the interior finite domain is sufficiently smooth. 

A similar approach has been successfully implemented by Jasiulek in 3D (|20|). 
Jasiulek's construction is based on Kerr-Schild coordinates and on a delicate control of 
characteristic speeds across the interface to the hyperboloidal domain. His numerical 
experiments with both finite difference and spectral methods in 3D show that his coordinates 
lead to more accurate evolutions than the spacelike matching of ( lfT7l[T8l ). In the following, 
we describe how transmitting layers based on Boyer-Lindquist coordinates are constructed. 

The Kerr metric in Boyer-Lindquist coordinates (t, r, 6, ip) reads 

/ 2Mr\ , o AaMr „ „ , , S , , 
5b l = - ( 1 - J dt 2 — sm 2 6 dtd<p+ — dr 2 

2 . 



+ S dd 2 + \ r ' z + a 2 + - ™ m " j sin 2 d dy? , (1) 

where X := r 2 + a 2 cosO 2 , and A := r 2 + a 2 — 2Mr. We denote the mass of the Kerr 
spacetime by M, its specific angular momentum by a. It is common to introduce the tortoise 
coordinate by 

dr* _ r 2 + a 2 

dr ~ A ' { ) 

A hyperboloidal (or more generally, transmitting) layer combines spatial compactifi- 
cation with a time transformation such that outgoing waves can be resolved near infinity. 
The spatial compactification maps the infinite domain E (-co, oo) to a finite domain 
p € [-S, S]. In this paper we use 

r* = £, with ft = l-(fc|:) 0(|p|-R), (3) 

where 9 denotes the step function and R the location of the interface. The function £1 vanishes 
at p = ±5 with a nonzero gradient. It is unity for — R < p < R and sufficiently smooth at 
the interfaces p = ±R. 

Such mappings from infinite to finite coordinate domains have been investigated as a 
potential solution to the outer boundary problem. Grosch and Orszag found that spatial 
compactification alone leads to resolution loss for hyperbolic problems ( 11271 ). However, 
outgoing waves can be resolved near infinity if a suitable time transformation is combined 
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with spatial compactification ( l23l ). Leaving the timelike Killing field invariant, the local 
time coordinate in the layer has the form 

T = t±h(r,). (4) 

The function h(r*) is called the height function and depends on the tortoise coordinate 
only. The time transformation is such that the asymptotic coordinate expression for the 
radially outgoing null direction is left invariant in the local coordinates of the layer (|28|). 
The Kerr metric in Boyer-Lindquist coordinates is asymptotically Schwarzschild in standard 
Schwarzschild coordinates. Therefore, the leading order form of in- and outgoing null rays 
can be written as t ± r*. We ignore the angular dependence because it falls off sufficiently 
fast ([24 1). As opposed to the characteristic formulation, only the asymptotic form of the 
null surfaces is important for the hyperboloidal transformation. This flexibility is the essential 
simplifying feature of hyperboloidal slicing that enables its straightforward application in Kerr 
spacetime. 

The invariance condition on the time function reads t ± r* = r ± p. Combined with the 
spatial compactification ^ and the form of the time transformation Q we get 

Hp) = ^-p- (5) 

As opposed to Ref. |24| where a hyperboloidal layer is attached in the positive direction 
of the tortoise coordinate only, here we attach hyperboloidal layers both in the positive and 
the negative directions thereby removing artificial boundaries completely from the numerical 
simulation (see also Ref. |25|). 



2.2. Numerical implementation 

We study the decay rates by solving the scalar wave equation 

□</> = 0, 

for a rescaled variable rip, using transmitting layers. The equations resulting from the formal 
transformations have been presented in the literature ( |[T8l 12410 . to which we refer the 
interested reader (the expressions <|3j and Q are inserted into the formal transformations). 

All our numerical simulations were performed with a modified version of the time- 
domain Teukolsky equation evolution code presented in Ref. 11 Oil . The code is an explicit, 
hyperbolic-PDE solver in (2+l)D that uses the Lax-Wendroff finite-difference evolution 
scheme. The numerical results have low truncation error due to high-order differencing in 
the angular direction and high-precision (quadruple or octal) floating-point operations. The 
code is parallelized via message passing interface using the standard domain-decomposition 
approach on the radial coordinate grid. More details on the code may be found in Ref. 1 10 1. 
The main modification of the code for this paper is the implementation of transmitting layers 
as discussed in the previous section. 

Typical grid resolutions used in this work are M/64 in the radial and 7r/48 in the angular 
direction. We set M = 1 and a = 0.995. The domain size in the compactified p-coordinate is 
[-100,100] with the transmitting layers beginning at R = ±28. The time-step is set to 1/128 
as dictated by the Courant condition for stability of the numerical evolution scheme. 

The initial data for our case studies are given by a generic (nonstationary) Gaussian 
wave-packet of width 4 centered at = 25. The angular dependence of the data is set by a 
pure, axisymmetric, spherical harmonic £'-mode. 
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3. Results 

We investigate the decay rates for all £' < 5 including the rates for the projected modes and 
rates for the full field to search for a possible appearance of splitting. We study intermediate 
decay rates for various observers by using the local decay rate n b s (t) defined as ([29]) 

dln\<f>(t,r ohs )\ 

n bs = -j-. 7 ■ (6) 

dint 

Note that the local decay rate is invariant under the transformations in the hyperboloidal layer 
because the timelike Killing field is left invariant in the layer. 



3.1. Far-field decay of projected modes 

The intermediate decay rates for projected modes with initial data £' < 5 measured by far-field 
observers ranging from null infinity to about 100AT are shown in Fig.[T] 

It is instructive to contrast Fig. [T] to observations in Schwarzschild spacetime. The 
observer dependence and the intermediate behavior of local decay rates are simple there. An 
initial pure spherical harmonic stays a pure spherical harmonic during evolution (£ — I') and 
the decay rates are given as Uj+ = — (2^ + 3) at timelike infinity and as n j? = —(£ + 2) at null 
infinity [1 13 j. During the intermediate decay there is a monotonous transition depending on 
the location of the observer f lfmPPHn . 

In Kerr spacetime, there is no geometric notion of a pure multipole because of the lack 
of spherical symmetry: the scalar spherical harmonics are not eigenfunctions of the Laplace 
operator, and this leads to coupling between multipoles. Therefore, even the question of the 
decay of pure initial £' modes may be considered artificial. Nevertheless, we observe a similar 
intermediate behavior in the projected modes for far field observers (r > lOOAf) in Fig.[T] 

The observer-dependence in these plots can be understood by considering the description 
of pointwise decay estimates suggested in Schwarzschild spacetime for small solutions 
( ||2T1 (221). For the dominant £' — £ = mode the estimate in standard Schwarzschild 
coordinates reads 

C 

(C2 + t + r*)(C 3 + t - 
where Cj are constants. This estimate captures the asymptotic behavior of the solution near 
null infinity as well as at finite distances from the black hole. It also provides an explicit 
description of the point- and slicing-dependence of the decay rates [ 12 1. The powers in the 
above formula are different for different values of £, but the qualitative structure of observer- 
dependence is the same. The far-field local-in-time decay rates suggested by Fig.[T]are listed 
in Tab. □ 
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Table 1. Far-field decay rates for all mode projections suggested by Fig. ^ We list both 
n j: and n i -|_ for each choice of £' ,£. The values agree with theoretical predictions: at null 
infinity nja- = -I' for t < I! - 2 and n.<? = -(£ + 2) for £ > l'\ at finite distances 
n i+ = - (I' + 1 + 1) for i < I' and n i+ =-(£' + £ + 3) for £>£' (see 0[lO)). However, 
some of the finite distance rates with £' = £, typed in boldface, may change asymptotically in 
time due to propagation of near-horizon splitting behavior omitted by Fig.[T](see Sec. |3.2fr . 
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Figure 1. Local decay rates for the evolution of even initial modes I 1 = 0, 2, 4 on the left 
panel and odd initial modes £' = 1, 3, 5 on the right panel, for a family of far-field observers 
ranging from null infinity to about lOOAf successively closer to the Kerr hole from top to 
bottom for each projection. The late-time values are consistent with the predictions in the 
published literature. No splitting is observed in the far field (r > lOOAf) for lOOOOAf. 



3.2. Near-field splitting in projected modes 

The picture of decay rates in certain projected modes presented in Fig.[T]and Tab.[T]is modified 
when the near horizon rates are taken into account. Racz and Toth report splitting in Ref. 1 19 [|] 
in the projected modes for certain cases with £' = £. We plot the near horizon decay rates for 
two cases, £' = I = 2 and f = I = 3 in Fig. [2] 

| We thank Gabor Zs Toth for pointing out that splitting occurs near the horizon for the projected modes. 
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Figure 2. Local decay rates for near-horizon observers with £' = i = 2 (left panel) 
and I 1 = I = 3 (right panel). The observers on the left panel are located at r*/M = 
20, 18, 16, 14, 12, 10 (yellow dashed), r*/M = 8, 6, 4, 2, 0, -4 (green solid) and r»/M = 
—8, —12, —20, —70, — co (blue dotted). The observers on the right panel are located at the 
same positions. The yellow dashed lines are consistent with Fig.^and Tab.^ The green solid 
vertical lines indicate the transition during which the field changes sign. The blue dotted lines 
are the near-horizon rates. The plot indicates that at asymptotically late times the near-horizon 
rates (blue) will dominate for far away observers (yellow) after a sign change (green). 



Decay rates along > 10M denoted by the yellow dashed lines are consistent with 
Fig.[T]and Tab.[T] The local rate curves bend down, however, in preparation for a sign change, 
which we see happening along the green solid lines in the range 6 [-4M, 8M]. The rate 
after the sign change approaches the same rate as the blue dotted lines, plotted for observers 
with 7% < — 8M. The green vertical lines in Fig. ^indicating sign change can be regarded as 
the boundary of splitting between the two finite distance decay rates. We see that the boundary 
moves farther out in time, but its motion is very slow. In fact, it seems to get slower with time. 

The figure suggests that asymptotically in time the decay of the blue dotted lines will 
dominate for all finite observers but at any given late time, there will be some far away 
observers that see the decay of the yellow dashed lines consistent with Fig. [Tj Note that we 
distinguish late time behavior for large t from asymptotic behavior for t — > oo. In that sense, 
splitting in the projected modes should be understood as an intermediate late time behavior, 
and not necessarily an asymptotic behavior. 

One can consider the following schematic model for the time-dependence of the field that 
leads to the switch in time. The field can be written as <j)(r, t) ~ A(r) t~ n + B(r) t~ n ~ l + 
C(r) t~ n ~ 2 + 0(t~ n ~ 3 ) (with possible additional logarithmic terms after the first two terms 
on the right hand side QUI ). The data suggest that A(r) x C(r) < and |C(r)| is large 
compared with |A(r)| for observers near the horizon, while |S(r)| is negligible. It would be 
interesting to see whether the sign difference and the relative strength of these terms can be 
computed using analytical methods, thereby explaining the near-horizon splitting. 

3.3. Far-field splitting in the full field 

A different type of splitting occurs in the decay rates of the full field as observed in Ref. 1 18 1. 
We plot in Fig. [3] the decay rates for the full field with f! = 4 (left panel) and = 5 (right 
panel). Consider the case I = 4. For a certain set of asymptotic observers the decay rates go 
below —5. For observers closer to the horizon, the rate curves up and approaches —5. 

To further emphasize the behavior of the full field, we present in Fig. |4]the local decay 
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Figure 3. Local decay rates for the evolution of £' = 4 (left panel) and I 1 = 5 (right panel) 
initial data. We observe splitting behavior in the full field, whereas no splitting occurs in 
the projected modes (compare FigJTJ. The splitting is explained by the relative strength of 
generated modes in Fig.|3]and Fig.|6j The lowest generated mode dominates along the yellow 
curves, whereas the higher mode is dominates along the Green ones. 
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Figure 4. The local decay rates for £' = A initial data as functions of inverse time M/t. The 
bold solid curve corresponds to an observer at r = 700Af , dashed curves for observers at 
r > 700Af, and thin solid curves to observers at r < 700M. For high M/t the distance 
of the observer increases from bottom to top. Note that the uppermost two thin solid curves 
would appear below the thin solid curves at the lower right corner were the figure extended 
to higher values of M/t. The dotted curves are extrapolations of the local decay rates using 
linear extrapolation for r < 700Af and quadratic extrapolations for r > 700Af. 



rate for the full field for initial data for £' — 4 as functions of inverse time for a family of 
observers. The curves separate into two types: those corresponding to near observers for 
whom the local decay rate at the end point of the numerical simulation increases as inverse 
time decreases, and those corresponding to far observers for whom it decreases. (Notice 
that the designations 'near' and 'far' here are different than in the preceding section.) The 
designation of near or far observers depends on how long the simulation runs. For our choice 
of ifinai = 10 4 M, the two types of curves separate at about r ~ 700A/. For near observers 
(r < 700Af) linear extrapolations of the local decay rate as a function of inverse time to 
M/t — > produce results very close to the expected asymptotic value of —5. More accurately, 
the observer at = 36 M measures an asymptotic decay rate of —5.017, the observer at 
= 65M measures an asymptotic decay rate of —5.028, the observer at = 112M 
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measures an asymptotic decay rate of —4.962, the observer at = 523A/ measures an 
asymptotic decay rate of —5.008, and that at 7% = 608M measures an asymptotic decay rate 
of —5.121. The observer at = 728AT measures an asymptotic decay rate of —5.017. Far 
observers have asymptotic values which increasingly deviate from the value of —5, as the 
simulation has not run sufficiently late to enter the asymptotic domain. However, there is 
enough curvature in the local decay rates of some far observers, such that the extrapolated 
values get close to —5. The observer at = 907Af measures an asymptotic decay rate 
of —5.058. More distant observers measure asymptotic rate of —5.495 for the observer at 
= 1,207M and —6.08 for the observer at = 1,806A/. However, notice that the 
latter two values vary significantly from the —5 value only because the simulation has not 
run far enough, such that the local decay rate do not include sufficient information about 
the asymptotic domain. Longer evolutions would make these observers near ones, and we 
predict that their asymptotic decay rate would be numerically close to —5. Notably, already 
with the given evolution time the extrapolations for the distant observers curve "up" in 
Fig. [4] in support of our prediction and explanation. Lastly, the observer at measures 
an asymptotic decay rate of —4.035. All these asymptotic values are numerically consistent 
with the expected integral values. 



^=0 ^=1 




Figure 5. Absolute values for the evolution measured by the observer along null infinity and at 
200Af for each projected mode. The lowest excited mode is the dominant mode. It is larger at 
all times by many orders of magnitude for I' =0,1 whereas for I 1 = 4, 5 the lowest excited 
mode dominates only at late times (see Fig. [6] for the ratio between the two lowest excited 
modes). The transition is observed as an intermediate splitting in the full field. 

Splitting of the full field is related to the relative strengths of modes generated due to 
mode coupling in Kerr spacetime. We plot in Fig. [5] the absolute values of the projected 




modes for initial data with £' = 0,1 on the top and £' = 4, 5 on the bottom panel. The values 
are measured by two observers for each mode: one at null infinity, the other at around 200A/. 

The dominant modes for £' = 0,1 are by many orders of magnitude larger than the up- 
modes generated by mode coupling (top panel of Fig.[5]l. This implies that the full field clearly 
decays with the rate of the dominant mode, and no splitting occurs. This clean separation of 
absolute values in projected modes is no more valid for initial data with £' = 4,5 as seen in the 
bottom panels of Fig. [5] The initial mode starts strongest, but decays fast leaving the generated 
lower modes behind. The null infinity decay rates of the lowest two modes are the same, 
implying that there is no transition of strength between them (compare also the lowest panel 
of Fig. [TJ. The finite distance rates of these modes, however, are different. Consequently, 
the lowest mode with slowest decay wins at late times, but at early times the decay of the 
faster decaying higher mode may still be valid. This leads to the splitting behavior at far away 
distances plotted in Fig. [3] 

Another way of demonstrating how splitting occurs in the full field is by looking at the 
ratio between the two lowest decaying modes plotted in Fig. [6] Because the rates of these 
modes along null infinity are the same, the decay along null infinity has no transition. Moving 
closer to the black hole, there are observers for which the ratio passes through unity. The 
passage through unity describes the transition of the decay rate from the rate of the higher 
mode to the lowest one. As indicated by Fig. [6] the finite distance observers will eventually 
observe the decay rate of the lowest mode. The far away the observer the longer it takes for 
the lowest generated mode to dominate. In that sense the splitting of the decay rates for the 
full field is only an intermediate behavior due to the relative strength of generated modes. 

4. Conclusions and Discussion 

We studied splitting of local decay rates observed in previous work ( lfl8l [191 l20l ). We 
discussed splitting in Kerr spacetime both in certain projected modes and in the full field. 
These two aspects of splitting seem unrelated to each other. Splitting in the projected modes 
appears near the black hole, and can be observed with standard codes. It is surprising that it 
has been discovered only recently. Splitting in the full field appears far away from the black 
hole and is due to the competition between the amplitudes of projected modes. 

Consider the example of £' = 4 initial data where full-field splitting appears (Fig. [3] left 
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panel). We observe that the initial mode decays the fastest. The generated lower modes 1 = 2 
and I = have the same decay rates at null infinity but different ones at finite distances. The 
I = mode decays slower at finite distances than the I = 2 mode. So, even though initially 
the generated I = 2 mode is stronger, eventually the I = mode dominates because of its 
slower decay (Fig. |5]>. The transition happens at different times for different observers, which 
appears as splitting in the decay rates of the full field. Asymptotically in time, the theoretical 
decay rates are valid. In that sense, the splitting of the full field is only an intermediate 
behavior. 

Splitting in certain projected modes with I' = £ is of a different nature. There we 
observe that for any given time there are three domains with three different local decay rates 
(see Figs. [T] and [2}: very close to the horizon, far away from the horizon, and near infinity. 
Figure [2] indicates that there is a transition in which the field changes sign and the decay rates 
approach the near-horizon rates. The boundary of this transition (denoted by vertical green 
lines in Fig. [2} moves very slowly towards infinity, which suggests that asymptotically in time 
the near-horizonrates will dominate at all finite distances. For all practical purposes, however, 
splitting is a real effect for these modes at any given (late) time. 

Note that the problem of the decay of prepared pure ^'-modes and its projected rates 
can be regarded rather artificial, and it is unlikely that splitting will have relevant physical 
consequences in the near future. Nevertheless it demonstrates a curious feature of mode 
coupling in Kerr spacetime. One would expect such behavior whenever there is mode 
coupling, also beyond the concrete example of Kerr spacetime. Therefore, we believe that 
this subject deserves an in-depth analytic study, for example by studying the branch cut in 
frequency domain (possibly similar to the detailed analysis in Schwarzschild spacetime of 
Ref. (3D). 

Finally, we mention that the transmitting layer technique presented in this paper should 
be useful for numerical time evolutions of dynamical fields in Kerr spacetime based on Boyer- 
Lindquist coordinates. Some interesting problems where this method might be used are 
superradiance (see, for example, recent work l32l ). and the instability of perturbations of 
extremal Kerr spacetimes [33]. 
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